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Abstract 

This paper presents a novel framework for track fitting which is usable in a wide range of experiments, independent of 
the specific event topology, detector setup, or magnetic field arrangement. This goal is achieved through a completely 
modular design. Fitting algorithms are implemented as interchangeable modules. At present, the framework contains 
a validated Kalman filter. Track parameterizations and the routines required to extrapolate the track parameters and 
their covariance matrices through the experiment are also implemented as interchangeable modules. Different track 
parameterizations and extrapolation routines can be used simultaneously for fitting of the same physical track. Rep- 
resentations of detector hits are the third modular ingredient to the framework. The hit dimensionality and orientation 
of planar tracking detectors are not restricted. Tracking information from detectors which do not measure the passage 
of particles in a fixed physical detector plane, e.g. drift chambers or TPCs, is used without any simplifications. The 
concept is implemented in a light-weight C++ library called GENFIT, which is available as free software. 
Keywords: track fitting, track reconstruction, Kalman filter, drift chamber, TPC 

1. Introduction 

Spectrometers in nuclear and particle physics have the purpose of identifying the 4-momenta and vertices of 
particles stemming from high-energy collisions and decays of particles or nuclei. In addition to calorimetric and other 
particle identification measurements, the 3-momenta and positions of charged particles are measured by tracking 
them in magnetic fields with the use of position sensitive detectors. Cluster finding procedures can be applied in some 
detectors to combine the responses of individual electronic channels in order to improve the accuracy of the position 
measurements. The position measurements will be called hits throughout this paper, regardless of whether they stem 
from a single detector channel or from a combination of several of them. Pattern recognition algorithms determine 
which hits contribute to the individual particle tracks. The hits identified at this stage to belong to one track then 
serve as the input for a fitting procedure, which determines the best estimates for the position and momentum of a 
particle at any point along its trajectory. A novel framework for this task of track fitting in complex detector systems 
is presented in this paper. It organizes the task of track fitting, i.e. the interplay between fitting algorithms, detector 
hits, and particles trajectories, with a minimal amount of interfaces. It is a toolkit which is independent of specific 
detector setups and magnetic field geometries and hence can be used for many particle physics experiments. 
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20 Tracking of particles is usually performed with a combination of different species of detectors. They can be cate- 

21 gorized according to the different geometrical information they deliver: 

22 1) detectors which measure the particle passage along one axis in a detector plane, e.g. silicon strip detectors or mul- 

23 tiwire proportional chambers; 

24 2) detectors which measure the two-dimensional penetration point of a particle through a plane, e.g. silicon pixel 

25 detectors; 

26 3) detectors which measure a drift time relative to a wire position, i.e. a surface of constant drift time around the wire 

27 through which the particle passed tangentially, e.g. drift chambers or "straw tubes"; 

28 4) detectors which measure three-dimensional space points on particle trajectories, like time projection chambers 

29 (TPC). 

so But also higher dimensional hits can occur: 

31 5) detector systems which measure two-dimensional position information in combination with two-dimensional di- 

32 rection information, including correlations between these parameters. Examples could be stations of several planes of 

33 detectors of categories 1 and 2, or electromagnetic calorimeters. 

34 For detectors which do not deliver tracking information in physical detector planes, e.g. those of categories 3 and 

35 4, the track fitting software of many experiments resorts to simplifications, which may be justified for a particular 
so application but prevent the usage of the same program for different experimental environments. Examples are the 
37 projection of TPC data onto planes defined by pad rows or the projection of the surfaces of constant drift time in drift 
sa chambers onto predefined planes, just leaving two lines with left-right ambiguities. This approach is problematic if 

39 the drift cells are not arranged in a planar configuration and if there is no preferred direction in which the detector is 

40 passed by the particles. Another common simplification is the treatment of two-dimensional hits (e.g. from silicon 

41 pixel detectors) as two independent one-dimensional measurements. 

42 In the framework presented here these problems have been overcome to make optimal use of the information from 

43 combinations of all types of tracking detector systems. All detector hits are defined in detector planes. For hits in 

44 detectors which do not have physical detector planes, so-called virtual detector planes are calculated dynamically for 

45 every extrapolation of a track to a hit. The dimensionality of detector hits is not restricted. One-dimensional hits 

46 constrain the track only along the coordinate axis in the detector plane which they measure. Two-dimensional hits 

47 are used in one fitting step to constrain the track in two dimensions in their detector planes. For hits in non-planar 

48 detectors (categories 3 and 4), the hit information (e.g. a surface of constant drift time) is converted into a position 

49 measurement in a plane perpendicular to the track, so that a fit is able to minimize the perpendicular distances between 
so the track and the position measurements. The information from hits with higher dimensionality, like those of category 

51 5, is used in four-dimensional hits, which contain all correlations between the parameters. 

52 Tracks of charged particles in magnetic fields are (usually) described by five parameters and a corresponding covari- 

53 ance matrix. The ability to extrapolate a track described by these parameters and their covariances, taking into account 

54 the effects of materials and magnetic fields, to different positions in the spectrometer is mandatory for track fitting. 
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The concept presented here provides a well defined interface for the invocation of external programs or libraries to 
perform these track extrapolations. It thus allows the straightforward use of established track following codes with 
their native geometry and magnetic field interfaces, such as GEANE [1|, which is nowadays distributed as part of 
CERN's Virtual Monte Carlo (VMC) package |2|. This is the most significant difference to other projects (e.g. Rec- 
Pack [3]), which offer more monolithic approaches to track fitting (e.g. defining their own geometry classes). The 
concept allows the simultaneous fitting of several representations of tracks to the same set of hits, i.e. to the same 
physical track. This flexibility is especially useful in the early phase of an experiment when different track parame- 
terizations and extrapolation approaches can be compared with each other, in order to identify the ones with optimal 
performance. But also the flexible coverage of different phase space regions with different track models, or the fitting 
of different mass hypotheses with the same track model can be desirable. The implementation of the concept has been 
realized in a software toolkit called GENFIT. It is written in C++ and is designed in a fully object oriented way. It 
has been developed in the framework of the PANDA experiment [4|, as part of the computing framework PANDAroot 
[5 1, but is now distributed as a stand-alone package 0. 

GENFIT contains a validated Kalman Filter. This algorithm is commonly used for track fitting in particle spec- 
trometers [7], since it performs much better than global minimization approaches in the presence of materials and 
inhomogeneous magnetic fields. The concept is however not limited to the use of the Kalman Filter. Other fitting 
algorithms, like Gaussian Sum Filters [8 1 or Deterministic Annealing Filters (5), can be implemented easily. 

Section [2] describes the concept of this new approach to track fitting in detail. Section [3] points out the key 
features of the implementation of GENFIT. Some examples of concrete track representations, on the dimensionalities 
of reconstruction hits and track representations, and the interplay between them follow in Sec. [4] Simulation studies 
which validate the Kalman filter implemented in GENFIT are presented in Sec. [5] 

2. Concept 

The basic functionalities which are required for any procedure of track fitting are the extrapolation of tracks to the 
positions of the hits in the detectors, and the calculation of the distances between hits and tracks, i.e. the residuals. The 
concept discussed here divides the problem of track fitting into three main entities which are separated from each other 
as much as possible and interact through well defined interfaces: 1) track fitting algorithms, 2) track representations, 
and 3) reconstruction hits. Figure [T] shows this structure. The following sections explain these objects in detail. 

2.1. Track Fitting Algorithms 

"Progressive" fitting algorithms like the extended Kalman filter (7] [10) are widely used for track fitting in high 
energy physics experiments. Although the track fitting concept discussed in this paper is not limited to the use of the 
Kalman filter, this algorithm shall serve as an example to illustrate which functionalities are generally required. 
The extended Kalman filter is an efficient recursive algorithm that finds the optimum estimate j& for the unknown true 
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Figure 1: General structure of objects for track fitting: Fitting Algorithm, Track Representation, and Reconstruction Hit. The arrows indicate the 
interactions between the objects, which are described in this chapter. POCA stands for point of closest approach. 



state vector % of a system from a series of noisy measurements, together with the corresponding covariance matrix C k 
of x k . The state vector contains the track parameters and the index k indicates that the state vector, and its covariance 
matrix are given at the detector plane of hit k. 

Before a recursion step, the state vector x K -\ an d covariance matrix C k -\ contain the information of all hits up to 
index k - 1 . In the prediction step the state vector and covariance matrix are extrapolated to the detector plane of 
hit k by the track following code. The predicted state vector is denoted by x k and the predicted covariance matrix by 
C k . This covariance matrix is the sum of the propagated track covariance matrix C k -\ (Gaussian error propagation by 
transformation with the Jacobian matrix of the propagation operation x k = /04-i))> an d a noise matrix which takes 
into account effects like multiple scattering and energy loss straggling. Then, the algorithm calculates the update for 
the state vector and the covariance matrix by taking into account the measurement m k : 

x k =~x k + K k ? k (1) 
C k = (I - K k H k )C k (2) 
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with the residual 



A = ™k - H k x k , (3) 

the weight of the residual (or Kalman gain) 

K k = C k H T k (H k C k H T k + V k y\ (4) 

and the covariance matrix V k of the measurement m k . / is the unit matrix of corresponding dimensionality. The 
projection matrix H k is a linear transformation from the coordinate system of the state vector x k , to the coordinate 
system of the position measurement m k of hit k, i.e. the detector plane of the hit. A discussion about dimensions of the 



as vectors and matrices in the above equations can be found in Sec. 4.2 together with concrete examples for the matrix 

87 H k . The elements of the covariance matrix C k shrink with the inclusion of more hits, thus reducing the impact of a 

as single hit on the value of the state vector. The \ 2 -contribution of hit kisx\ - ^(Vk ~ H k C k H^Y l r k with the filtered 

89 residual r k — m k - H k x k . It adds dim(m^) degrees of freedom to the total x 1 - 

90 After the Kalman steps have been performed on all hits of the track, the track can still be biased due to wrong starting 

91 values x®. This bias can be reduced by the repeated application of the procedure in the opposite order of hits, using the 

92 previous fit result as starting values for the track parameters. Before the fit is repeated, the elements of the covariance 

93 matrix have to be multiplied with a large factor ((9(1000)) in order not to include the same information in the track 

94 several times. 

95 As can be seen in Fig.[T]the fitting algorithm operates on entities called reconstruction hits and track representations, 

96 which are detailed in the following. 

97 2.2. Track Representations 

98 A particle track is described by a set of track parameters and a corresponding covariance matrix, which are defined 

99 at a given position along the track. Often, the track parameters are e.g. given at a particular z-position. In the concept 

100 presented here, track parameters are always defined in reference planes. 

101 In order to use a track model in a track fitter, one needs to be able to extrapolate the track parameters to different 

102 places in the spectrometer. The combination of the track parameterization and the track extrapolation functionality 

103 will be called a track representation. A track representation holds the data of the state vector x k , and the covariance 

104 matrix C k of a track, as well as the reference plane at which these are defined. Also it provides a well defined interface 

105 for the invocation of the external routines needed to extrapolate the parameters to different positions. As can be seen 

106 in Fig.[T] there are three track extrapolation functions which are needed for each track representation: Extrapolation 

107 to a plane, extrapolation to the point of closest approach (POCA) to a point, and extrapolation to the point of closest 

108 approach to a line. Fitting algorithms access the track parameters and extrapolation functions in a common way via 

109 the track representation interface without knowledge of the specific form of the track parameterization or the way the 
no extrapolations are carried out. 
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in Different track representations can be used in parallel. It is possible to fit the same track, i.e. the same set of hits, with 

112 different track representations simultaneously. There are several reasons why this is desirable: For low momentum 

113 particles the fitting of different mass hypotheses with the same track representation can give a clue to the particle 

114 identity via the^ 2 of the fits, because the different energy loss for different particle masses at a given momentum leads 

115 to different extrapolations. Fitting of the same track with different parameterizations and extrapolation tools can be 
us advantageous as well. In the early phase of an experiment one can compare different track representations to identify 
117 the ones which perform best, or there could be regions in phase space in experiments where it might not be clear 
us beforehand which track representation will give the best results. Then one can just fit several of them simultaneously 

119 and retain the best result. 

120 2.3. Reconstruction Hits 

121 The object which represents a position measurement from a detector used in a track fit is called a reconstruction hit. 

122 It contains the vector of the raw measurement coordinates and its corresponding covariance matrix. As discussed in 

123 the introduction, the nature of this raw hit information can be quite diverse. It can e.g. be a direct position measurement 

124 or a drift time. As can be seen in Fig.[T] a reconstruction hit provides its detector plane, the measurement coordinates 

125 mj. in the detector plane coordinate system, the covariance matrix in the detector plane coordinate system, and 

126 the projection matrix to the fitting algorithm. For detectors, which measure positions in a physical detector plane 

127 (categories 1 and 2 of Sec.[T]i, the detector plane is identical with the physical plane. 

128 For non-planar detectors like wire-based drift chambers or TPCs (categories 3 and 4 of Sec. [TJ, no such physical 

129 detector planes are defined. Instead, the concept of virtual detector planes is introduced. For space-point detectors, 

130 the track fit has to minimize the perpendicular distances of the track to the hits. Therefore, the virtual detector plane 

131 for each hit must contain the hit position and the point of closest approach of the track to the hit point. Then the 

132 residual vector which points from the hit point to the point of closest approach will be perpendicular to the track. This 

133 geometry is illustrated in Fig. [5] The orientation of the spanning vectors it and v is chosen arbitrarily in the plane. For 

134 wire-based drift detectors the virtual detector plane contains the point of closest approach of the track to the wire, and 

135 is oriented to contain the whole wire. The spanning vectors are chosen to lie perpendicular (it) and along (v) the wire. 

136 This geometry is shown in Fig. [3] The wire position and drift time are then measurements of u (the v coordinate could 

137 be measured via double-sided readout with charge sharing or time of propagation). In both cases, the orientation of 

138 the plane will directly depend on the track parameters. The consequence is that virtual detector planes have to be 

139 calculated each time a hit is to be used in a fitting step. The reconstruction hit uses the corresponding extrapolation 

140 function of the given track representation to find the point of closest approach as indicated in Fig.[T] 

141 Different kinds of reconstruction hits are accessed via a common interface. When the fitting algorithm obtains the 
u2 detector plane from a reconstruction hit, it does not know whether it will receive a physical or a virtual detector plane. 

143 This distinction is fully handled inside the reconstruction hit. 

144 After the detector plane is defined, the reconstruction hit can provide the measurement coordinate vector m^, and the 
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145 hit covariance matrix For non-planar detectors, these quantities are results of coordinate transformations into the 

146 virtual detector plane (hence the difference between the raw hit coordinates/covariance and the vector ink and matrix 

147 Vj in Fig. [TJ. The three-dimensional hit vector and the 3x3 covariance matrix of a space-point hit are transformed 

148 into a two-dimensional vector in the detector plane and a 2 x 2 covariance matrix. Even if the errors of the space 

149 point were uncorrelated, the matrix V* will in general contain a correlation, which is taken into account in the fit. For 

150 wire-based drift chambers, the drift time information is converted to a position information in the calculation of ihk 

151 and Vk- 

152 The projection matrix Hk transforms the state vector from the given track parameterization into the coordinate system 

153 of the hit. In order to determine this matrix, the concrete coordinate systems of the track representation and the 

154 reconstruction hit must be known. Since there will be typically more different types of reconstruction hits than track 

155 representations, the projection matrix is determined in the reconstruction hit object. The matrix Hk provides the only 

156 link between a given track parameterization and the different hit coordinate systems. If a fit is performed with several 

157 track representations, the same reconstruction hit will provide a different matrix Hk for each track representation. 

158 3. Implementation - GENFIT 

159 The software package which implements the concept presented in this paper is called GENFIT J6). It is completely 

160 written in C++ and makes extensive use of object oriented design. It uses the C++ standard template library [11] and 

161 the ROOT data analysis framework [ 12|. 

162 Figure [4] presents the general class structure of GENFIT. The classes representing the fitting algorithms operate on 

163 instances of the class GFTracli^] A GFTrack object contains a std: : vector<GFAbsRecoHit*> and a 



164 std : : vector<GFAbsTrackRep*>. The reconstruction hits and track representations of Sees. 2.2 and 2.3 are realized 

165 as polymorphic classes. The class GFAbsRecoHit is the interface class to the reconstruction hits, and GFAbsTrackRep 

166 is the interface class to the track representations. 

167 The reconstruction hit objects are created from the position information acquired in the detectors. The pattern recog- 

168 nition algorithms, which precede the use of GENFIT, determine which of these detector hits belong to a certain track. 

169 They deliver an instance of the class GFTrackCand, which holds a list of indices which identify the hits belonging 

170 to the track. A mechanism called GFRecoHitFactory has been implemented to load the reconstruction hits into the 

171 GFTrack object. 

172 3.1. Track Representations 

173 In order to use a particular track parameterization for track fitting in GENFIT, one needs code which can ex- 

174 trapolate such track parameters, taking into account material effects on the track parameters and their covariance 

175 matrix. In order to interface the track model to GENFIT, one implements a C++ class which inherits from the 



2 class names or other code fragments are set in typewriter font. 
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178 



abstract base class GFAbsTrackRep and provides an implementation for the virtual methods extrapolate ( . . . ), 
extrapolateToPoint ( . . .), and extrapolateToLine( . . .). Section 4.1 presents examples of concrete track 
representations. 



179 3.2. Reconstruction Hits 

iao The fitting algorithms interact with the reconstruction hits via the abstract base class GFAbsRecoHit. The recon- 

181 struction hits do, however, not inherit directly from this class, but from the intermediate interface class 

182 GFRecoHitlf c<Policy>. This is illustrated in Fig. [5] For more information about the policy design pattern, please 

183 see ITT31 . There are (currently) three geometrical categories of reconstruction hits: Hits in planar detectors, space- 

184 point hits, and hits in wire-based drift chambers, which deliver their wire position and a drift time. This catego- 

185 rization is expressed in the code by the three different policy classes GFPlanarPolicy, GFSpacepointPolicy, and 

186 GFWirePolicy. These policy classes all implement functions for calculating or delivering the detector plane, the 

187 hit coordinates in the detector plane, and the hit covariance matrix in the detector plane. They are used to unify the 

188 geometrical properties of reconstruction hits to avoid any code duplication in the implementation of similar recon- 

189 struction hits. The latter two policies use the corresponding track representations to calculate the virtual detector 



planes, as detailed in Sec. 2.3 



191 As described in Sec. 2. 1 the fitting algorithm needs a matrix which is a linear transformation from the vector space 

192 of track parameters to the coordinate system defined by the detector plane. The virtual method 

193 GFAbsRecoHit : :getHMatrix( . . . ) is overridden in the implementations of the concrete reconstruction hits. In 

194 order to provide the correct matrix, the reconstruction hit determines the concrete typ^jof the track representation it 

195 is asked to interact with in this particular fitting step. This type checking is represented by the lower arrow in Fig.[T] 

196 It is the only place in GENFIT where a direct type compatibility of tracks and hits is checked. A maximal modularity 

197 of the system is achieved through this mechanism. If one adds an additional track representation, it is quite obvious 

198 that one has to provide new coordinate transformations from this new parameter space into the coordinate systems in 

199 which the hits are defined. 



200 4. Examples 

201 4.1. Concrete Track Representations 

202 A concrete interface to an external track propagation package which has been realized with GENFIT is the track 

203 representation called GeaneTrackRep2. It is based on the FORTRAN code GEANE. The detector geometry is in- 

204 eluded via the TGeo classes of ROOT lfT2l and the magnetic field maps are accessed via a simple interface class called 

205 GFAbsBField. State vectors for this track representation are defined as = (q/\p\,du/dw,dv/dw,u,v) T , where the 



3 by performing a C++ dynamic_cast on the base class pointer GFAbsTrackRep*. 
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206 detector plane is spanned by the vectors u and v (normal vector w = it x q denotes the particle charge and p is the 

207 particle momentum. The quantitative tests of GENFIT in Sec.[5]are carried out with this track representation. 

208 Another track representation included in the GENFIT distribution is called RKTrackRep. It was adopted from the 

209 COMPASS experiment |T4l and uses a Runge-Kutta solver to follow particles through magnetic fields. It has the 

210 same state vector definition as GeaneTrackRep2. It also uses the TGeo classes for the geometry interface. 

211 4.2. Interplay between Track Representations and Reconstruction Hits 

212 The classes which represent the fitting algorithms just carry out their linear algebra without knowing about the 

213 dimensions of the state vectors j4 and the measurement vectors fhk- The matrix Hk is provided by the reconstruction 
2u hit class to transform state vectors and covariance matrices of a specific parameterization into the measurement vector 

215 coordinate system. This projection matrix ensures that the dimensionalities of the vectors and matrices in the fitting 

216 algorithm are compatible with each other. The following examples shall illustrate this: 

217 1. A four-dimensional track model can be used for tracking without magnetic fields. The state vector is defined 

218 as Xk = (u, v, du/dw, dv/dw) T for a straight line where u and v span the detector plane, and w = u x v is the 

219 normal vector. A strip detector shall measure the u coordinate. Then the measurement vector of equation 

220 ink, is a scalar. The projection matrix is defined as Hk = (1,0, 0, 0), so that Hk ■ Xk is one-dimensional, just as 

221 the residual r^. The Kalman gain is a 4 x 1 matrix, and the x 1 -increment is correctly calculated for one degree 

222 of freedom, in the sense that f*k and (V* - HkCkH^) are scalars. 

2. A pixel detector is used in combination with a five-dimensional trajectory model for charged particle tracking 
in magnetic fields. The detector measures the coordinates u and v in the detector plane, and the state vector is 
■4 = (q/\P\> du/dw, dv/dw, u, v) T . The 2x5 projection matrix is then: 



H k = 



1 N 
1 



All matrices and vectors automatically appear with correct dimensions: ihk and f*k are 2-vectors, Vk is a 2 x 2 
matrix, the Kalman gain is a 5 x 2 matrix, and^ 2 is a scalar which is calculated from two degrees of freedom 
(rk is a 2-vector, and (Vk - HkCkH^) is a 2 X 2 matrix). 

If the next hit in the same track only measures one coordinate (e.g. u in the detector plane coordinate system of 
the next hit) , fhk will be scalar, Hk will be of dimension 1x5, and there will be only one degree of freedom 
added to the overall^- 2 . 

3. A TPC delivers space-point hits. The track model is the same as in example 2. The TPC measures three spatial 
coordinates but this information is transformed into a two-dimensional hit in the virtual detector plane, which 
is perpendicular to the track. This two-dimensional hit is treated identically to example 2. This is the desired 
behaviour, since measurements or errors along the flight direction do not contribute to the track fit. 
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233 5. Simulation Studies 



234 The statistical and numerical correctness of a Kalman filter fit depends on the following items: 1) The mathematics 

235 of the Kalman filter have to be implemented correctly. 2) The projections of the covariance matrices of the hits onto 

236 the (virtual) detector planes have to be correct. 3) The propagation of the track parameters and the covariance matrix 

237 are done correctly. For the covariance matrix this means the correct estimation of the Jacobian matrices needed for 

238 the Gaussian error propagation. 4) The effects of traversed materials must be taken into account correctly: the state 

239 vector has to be modified (momentum loss) and the entries of the covariance matrix need to be increased by the 

240 addition of noise matrices (e.g. due to multiple scattering) [7|. Since the track representations are external modules, 

241 the Kalman filter and the reconstruction hit implementation in GENFIT are tested with a simplified setup, where the 

242 particles traverse a vacuum. This way, the effects number 1 to 3 are tested, while the effect number 4 is decoupled 

243 and not tested here. The setup contains a homogeneous magnetic field, since possible problems arising from field 

244 inhomogeneities would only point to problems in the external track representation module and not in the GENFIT core 

245 classes. Instead of detector responses with full digitization simulations, which result in unknown detector resolutions, 

246 known measurement errors are used. 

247 The track representation GeaneTrackRep2 is used for these tests. The program samples 30 space points on the 

248 trajectory at distances of 1 cm, which are smeared with Gaussian distributions of known widths. Like in a TPC, 

249 the x- and y-measurements are assumed to have equal and better resolutions than the z-coordinate measurements 

250 (cr v = cr y = 1/2 • cr z ). These smeared points are used in the fit as reconstruction hits based on GFSpacepointPolicy 

251 similar to TPC measurements (see Pig . |5j . In front of the first hit, a reference plane is defined in which the fitted track 

252 parameters are compared to their true values to obtain residual and pull distributions^] If the fit is able to correctly 

253 determine the track parameters and their errors, the pull distributions will be Gaussians of width cr = 1 and of mean 

254 value 0. Figure [6] shows the five pull distributions for the track parameters, which fulfill these criteria within the 

255 corresponding errors, proving that the non-uniform errors of the hits are taken into account correctly. 

256 Another test is carried out with a slightly different detector geometry. Hits from 15 crossed planes of strip detectors 

257 are fitted together with 15 space-point hits. The strip hits each contribute one degree of freedom, the space-point hits 

258 each contribute two degrees of freedom (they only constrain the track in a plane perpendicular to the track), and the 

259 track parameters subtract five degrees of freedom (15 + 2 • 15 — 5 = 40). The ^-probabilities for these fits are shown 

260 in Fig. [7] If the number of degrees of freedom is taken into account correctly, this distribution is expected to be fiat. 

261 A;f 2 -test against a uniform distribution results in a^ 2 /n.d.f. = 87.1/99, close to unity, as expected. 

TM 

262 The execution time per track is 14 ms on one core of an AMD Phenom II X4 940 CPU for 30 space-point hits 

263 with one forward and one backward fitting pass of all hits. Of this time, a fraction of 91% is spent in the external 

264 extrapolation routines of GeaneTrackRep2, as determined with Valgrind lfT31 . The GENFIT core classes have not 

4 the pull of a variable x is defined as (xfl t - x Uue )/o- x . 
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266 



yet been optimized for execution time, but the above result shows that optimizations would be most rewarding in the 
track extrapolation routines. 



267 6. Conclusions and Outlook 

268 A novel framework for track fitting in particle physics experiments has been presented in this paper. Its implemen- 

269 tation is a C++ library called GENFIT, which is available freely. Its modular design consists of three major building 

270 blocks: Fitting algorithms, track representations, and reconstruction hits. GENFIT contains a validated Kalman filter. 

271 A standard Kalman smoother is planned to be implemented in the future, as well as other fitting algorithms. The 

272 possibility of the application of GENFIT to pattern recognition tasks seems promising and will be investigated. 

273 The generic design of the track representation interface enables the user to use any external track following code 

274 with GENFIT. The framework allows simultaneous fits of the same particle track with different track representations. 

275 Possible applications of this feature are the fitting of different mass hypotheses with the same track model, or the test 

276 and validation of different track parameterizations and track following codes. Also the coverage of different regions of 

277 phase space with specialized track representations is an important feature in many experiments. At present, GENFIT 

278 contains two track representations which provide interfaces to the track following code GEANE and a Runge-Kutta 

279 based track extrapolation code adopted from the COMPASS experiment. New track representations which allow the 

280 use of other track following codes can be implemented in a straightforward way. The interfaces to the detector geom- 

281 etry and the magnetic field maps can be chosen freely and are all encapsulated in the track representation class. 

282 The geometrical properties of reconstruction hits are not restricted in this framework. The dimensionality of hits is not 

283 fixed to particular values, and the orientation of detector planes can be chosen freely. Hits from detectors which do not 

284 measure the passage of particles in predefined planes, such as drift chambers or TPCs, are handled in the concept of 

285 virtual detector planes. This leads to a direct minimization of the perpendicular distances between the particle tracks 

286 and the position measurements from the detectors, i.e. the surfaces of constant drift time or the space points measured 

287 in a TPC. 

288 GENFIT provides an easy-to-use toolkit for track fitting to the community of nuclear and particle physics. It is used 

289 in the PANDA computing framework. Applications in other experiments are being considered (e.g. Belle II). 
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Figure 2: Virtual detector plane (spanning vectors it and for a space-point hit. 
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Figure 3: Virtual detector plane (spanning vectors it and \f) for a wire-based drift detector. 
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Figure 4: Class structure of GENFIT. The detailed inheritance structure of reconstruction hits is shown in Fig. [5] 
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Figure 5: Inheritance structure of reconstruction hits in GENFIT. 
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pull of variable 



Figure 6: Pull distributions for the five track parameters of GeaneTrackRep2. The pull of a variable x is defined as (xg t - x Uw )I<t x . 
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Figure 7: ^-probability distribution for a series of track fits through 15 planes of strip detectors and 15 space points. 



18 



